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1. Lattice setup 

We report about an ongoing project concerned with the computation of the spectrum of mesons 
with at least one strange or charm valence quark using twisted mass lattice QCD with 2+1+1 
dynamical quark flavors. Here we present first results for masses of ground state and excited D and 
Dg mesons. 

Currently we have performed computations on f» 400 gauge field configurations from a single 
ensemble generated by the European Twisted Mass Collaboration (ETMC). This ensemble is char- 
acterized by 24^ X 48 lattice sites, a lattice spacing a ^ 0.086 fm and a pion mass nijc ~ 324 MeV. 
For detailed information regarding these gauge field configurations we refer to [p. 

For both the valence strange and charm quarks we use degenerate twisted mass doublets, i.e. 
a different discretization as for the corresponding sea quarks. We do this, to avoid mixing of 
strange and charm quarks, which inevitably takes place in a unitary setup, and which is particularly 
problematic for hadrons containing charm quarks [^. The degenerate valence doublets allow two 
realizations for strange as well as for charm quarks, either with a twisted mass term +i7i.s,c75 or 
—iHsxYs- While for a qq meson creation operator the sign combinations (+, — ) and (— , +) for the 
quark q and the antiquark q are related by symmetry, (+,+) and (—,—), which are also related, 
yield at finite lattice spacing results slightly different compared to (+,— ) and (—,+), due to dif- 
ferent discretization errors. Using (+,— ) = (— ,+) we have tuned the valence strange and charm 
quark masses to reproduce the physical values of 2mj. — m\ and mo, two quantities, which only 
weakly depend on the light u/d quark mass. 

We have constructed suitable D and Ds meson creation operators using the maximum of sixteen 
independent y combinations as well as lattice versions of spherical harmonics. Moreover, APE 
and Gaussian smearing is used, to optimize the overlap of the trial states to the low lying energy 
eigenstates of interest. More details regarding the construction of meson creation operators in 
twisted mass lattice QCD can be found in 

For the computation of the corresponding correlation matrices we use spin diluted timeslice 
sources in combination with the "one-end trick". Meson masses are then determined from plateaux 
values of corresponding effective masses, which we obtain by "diagonahzing" the correlation ma- 
trices by solving generalized eigenvalue problems. 

2. Results 

In Figure |I| we show D and meson masses for quantum numbers 7^ = 0^,0+, 1^, 1+. Blue 
data points represent experimental results, while red and green data points correspond to our lattice 
results obtained with sign combinations (+,— ) = (— ,+) and (+,+) = (—,—), respectively. The 
differences between the "red and green lattice results" are rather small, indicating that discretization 
errors at the currently used lattice spacing are <2%. 

While in most cases lattice and experimental results agree within statistical errors, we observe 
significant discrepancies for D*,, and D^i. Similar findings have been reported by other lattice 
studies and by model calculations (cf. e.g. [Q, |5|, ^). These discrepancies support the assumption 
that D*Q and D^i are not ordinary quark antiquark states, but might e.g. be mesonic D-K molecules 
or diquark-antidiquark pairs (tetraquarks). A corresponding investigation of possibly present four- 
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Figure 1: Low lying D meson (left) and Dg meson (right) spectrum. 

quark components in positive parity mesons using the same lattice setup is ongoing [^]. A study of 
the broad Dq and Di mesons as resonances using lattice QCD and a composition of usual quark- 
antiquark states with D + k and D* + K can be found in 

In the near future we plan to extend our computations to states with spin 7 > 1 and to the 
charmonium spectrum. We will also use ensembles at even smaller lattice spacings to study the 
continuum limit and we will consider several different pion masses to investigate the light u/d 
quark mass dependence. 

Acknowledgments 

We thank Christian Wiese for discussions. M.K. and M.W acknowledge support by the 
Emmy Noether Programme of the DFG (German Research Foundation), grant WA 3000/1-1 and 
by Helmholtz Graduate School HGS-HIRe for FAIR. This work was supported in part by the 
Helmholtz International Center for FAIR within the framework of the LOEWE program launched 
by the State of Hesse. 



References 

[1] R. Baron et al. [ETM Collaboration], JHEP 1006, 1 1 1 (2010) [arXiv: 1004.5284 [hep-lat]]. 

[2] R. Baron et al. [ETM Collaboration], Comput. Phys. Commun. 182, 299 (201 1) [arXiv: 1005.2042 
[hep-lat]]. 

[3] K. Jansen et al. [ETM Collaboration], JHEP 0812, 058 (2008) [arXiv:0810.1843 [hep-lat]]. 

[4] D. Mohler and R. M. Woloshyn, Phys. Rev. D 84, 054505 (201 1) [arXiv:l 103.5506 [hep-lat]]. 

[5] G. Moir and M. Peardon, talk given at "Xth Quark Confinement and the Hadron Spectrum," October 
8-12, 2012, Munich, Germany. 

[6] D. Ebert, R. N. Faustov and V. O. Galkin, Eur. Phys. J. C 66, 197 (2010) [arXiv:0910.5612 [hep-ph]]. 

[7] J. O. Daldrop, C. Alexandrou, M. D. Brida, M. Gravina, L. Scorzato, C. Urbach and M. Wagner 
[ETM Collaboration], PoS LATTICE 2012, 161 (2012) [arXiv:121 1.5002 [hep-lat]]. 

[8] D. Mohler, S. Prelovsek and R. M. Woloshyn, arXiv: 1208.4059 [hep-lat]. 



3 



